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ABSTRACT 

We calculate the cross-correlation function (CCF) between damped Ly-a systems 
(DLAs) and Lyman break galaxies (LBGs) using cosmological hydrodynamic simula- 
tions at z = 3. We compute the CCF with two different methods. First, we assume 
that there is one DLA in each dark matter halo if its DLA cross section is non-zero. In 
our second approach we weight the pair-count by the DLA cross section of each halo, 
yielding a cross-section-weighted CCF. We also compute the angular CCF for direct 
comparison with observations. Finally, we calculate the auto-correlation functions of 
LBGs and DLAs, and their bias against the dark matter distribution. For these differ- 
ent approaches, we consistently find that t here is good agreem ent b etween our simu- 
lations and observational measurements bv lCooke et al.l ()2006af ) and lAdelberger et alj 
([2005). Our results thus confirm that the spatial distribution of LBGs can be well 
described within the framework of the concordance ACDM model, and support the 
argument that the distribution of DLAs is strongly correlated with that of LBGs. 

Key words: cosmology: theory — stars: formation — galaxies: evolution - galaxies: 
formation - methods: numerical 



1 INTRODUCTION 

According to the cold dark matter (CDM) model of struc- 
ture formation, the spatial distribution of galaxies can be 
understood as a result of the gravitational instability of den- 
sity fluctuations in the CDM, and the dark matter halo 
mass function can be well described by analytic models 
l|Sheth fc Tormen|[l999h . More precisely, hierarchical CDM 
models predict that the massive galaxies at high redshift 
(hereafter high-2:) are clustered together in high-density re- 
gions, while low-mass galaxies tend to be more evenly spaced 
(Sardeen et al.lll986l : lKaisej|l984 ). Under the assumption 
that galaxies are produced from primordial density fluctua- 
tions owing to gravitational instability, one can estimate the 
average mass of galaxy host ha l os ba sed on clustering data. 
For example, lAdelberger et all (|2003h estimated the typical 
halo mass of LBGs at 2; ~ 3 to be Mhaio — 10^^ Mq from 
observations of their auto-correlation function (ACF). 

Damped Lyman-Q systems (DLAs), defined as quasar 
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absorption sy stems with column density of A'hi > 2 x 10 
atoms cm~^ l|Wolfe et al.lll986l ). probe the Hi gas associ- 
ated with high-z galaxies. Since stars are hardly formed in 
warm ionized gas and are tightly correlated with cold neu- 
tral clouds, the amount of Hi gas is very important, be - 
ing the precursor of molecular clouds (IWolfe et al.ll2003l ). 
DLAs dominate the Hi content of the Universe at 2: ~ 3 
and contain a sufficient amount of Hi gas mass to ac- 
count for a large fraction of th e present-day stellar mass 
IStorrie-Lombardi fc Woi^l2000D . The gas kinematics and 
chemical abundances of DLAs can be measured and are doc- 
umented in detail. However, the masses of DLA host halos 
(hereafter DLA halos) remain poorly constrained, because 
only about a quarter of all z > 3 quasars exhibit DLA ab- 
sorption, and the scattered distribution of DLAs in quasar 
sight lines precludes the use of DLAs as tracers of dark mat- 
ter halo mass. 

Alternatively, the mass of DLA halos can be probed 
by the cross-correlation between DLAs and a galaxy pop- 
ulation whose clustering and halo mass are well under- 
stood. iCooke et al.. (, 2006a . b) used 211 LBG spectra and 
11 DLAs to measure the three dime nsional (3-D) LBG 
ACF and DLA-LBG CCF fsee also iBouche et al.l l2005l : 
iBouche fc Lowenthalll2004l : lGawiser et al.ll200l[) . Theh anaL 
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ysis started by counting the number of LBGs in 3-D cylindri- 
c al bins centered on each of 11 DLAs, following the method 
of lAdelberger et al.l (|2003l ). They detected a statistically sig- 
nificant result of DLA-LBG CCF, and estimated an average 
DLA halo mass of (A/dla) ~ 10^^'^ Mq, assuming a single 
galaxy per halo. 

On the theoretical side, iNagamine et al.l (|2007l ) calcu- 
lated the average DLA halo mass using a series of cosmo- 
logical hydrodynamic simulations with different box sizes, 
resolution and feedback strengths. They found a mean DLA 
halo m ass of (Mdla) = 10 ^^'^ Mp, which is comparable 
to what Cooke et al. feOOeal lbh obtained. More recent work 
by fPontzcn ct al. ( 200^) also showed that the DLA cross- 
section is predominantly provided by intermediate mass ha- 
los, 10^ < Mvir/M0 < lO". These results motivate us 
to further examine the distribution of DLAs relative to 
that of LBGs. In this paper, we compute the DLA-LBG 
CCF in cosmological SPH simu lations, using the sample of 
LBGs and DLAs obtained by iNagamine et al.l ()2004al lbh. 
We compare our result s wit h the observational results by 
lAdelberger et al.l (|2005l ) and lCooke eral] l|2006al lbh. 

Our paper is organized as follows. In Section (2] we 
briefly describe the features of our cosmological SPH sim- 
ulations used in this paper. In Section [3] and Section HI 
we describe and report the methodology, binning method, 
and the results for 'unweighted' and 'weighted' DLA-LBG 
CCF, respectively. We then discuss the projected angular 
CCF for the direct comparison with observational result by 
ICooke et al.l l|2006al lbl) in Sectiong] The ACFs of LBG-LBG 
and DLA-DLA are discussed in Subsections 16.11 and 16.21 
while the bias results are reported in Section [T] Finally, we 
discuss the implications of our work in Section [H] 



2 SIMULATIONS 

In this paper, we use two different cosmological 
smoothed particle hydrod ynamics (SPH) simulations 
l|Springel fc Hcrnguist 2003b') performed with the GADGET- 
2 code (fSprin gel 2005). The simulation parameters of the 
two runs (named D5 and G5) are summa rized in Table [1] 
The sa me set of runs has been used by INagamine et al.l 
l|2004al lbl. 120071 ) to study the global properties of DLAs, 
such as the DLA cross section, incidence rate, and Hi 
column density distribution functions. 

The code we use is characterized by four main features. 
First, it uses an entropy-conserving formulation of SPH 
l|Springel fc Herrigu ist 2002), which explicitly conserves en- 
tropy of the gas where appropriate. Second, highly overdense 
gas particles are treated with a sub-resolution model for th e 
interstellar medium (ISM) jSpringel fc Hernguis^ l2003al '). 
The dense ISM is here assumed to be made of a two-phase 
fluid consisting of cold clouds in pressure equilibrium with 
a hot ambient phase. Cold clouds grow by radiative cooling, 
and form the reservoir of baryons for star formation. Once 
star formation occurs, the resulting supernovae (SNe) de- 
posit energy into the ISM, heating the hot gas environment, 
evaporating cold clouds, and transferring cold gas back into 
the ambient phase. This established a self-regulation cycle 
for star formation in the ISM, Additionally, the simulation 
keeps track of metal abundance and the dynamical trans- 



Run 


i'box 


iVp 






e 


D5 

G5 


33.75 
100.0 


2 X 324^ 
2 X 324^ 


8.15 X 10^ 
2.12 X 10^ 


1.26 X 10^ 
3.26 X 10® 


4.17 
12.3 



Table 1. Simulations employed in this study. A'^p is the initial 
particle number of gas and dark matter particles (hence x2). 
™DM S'lid rrigas are the masses of dark matter and gas particles 
in units of h~^MQ, respectively, e is the comoving gravitational 
softening length in units of kpc, which is a measure of spatial 
resolution. All runs adopt a 'strong' galactic wind feedback model. 



port of metals. Metals are produced by stars and returned 
into the gas by SNe. 

Third, a model for galactic winds is included to study 
the effects of outflows on DLAs, galaxies, and the intergalac- 
tic medium (IGM). In this model, gas particles are driven 
out of the dense star-forming medium by assi gning an ex- 
tra mo mentum in random directions t Springel fc Hernguistl 
l2003al ). It is assumed that the wind mass- loss rate is pro- 
portional to the star formation rate, and the wind takes a 
fixed fraction of the SN energy. For the D5 and G5 runs, 
a strong wind speed of 484kms~^ is adopted. Fourth, the 
code includes radiative cooli ng and heating with a un iform 
UV b ackground of a modified Haardt fc Madaul l|l996l 'l spec- 
trum l|Dave et al.ll 19991 : iKatz et al.lll996l ). where reionisation 
takes place at 2; ~ 6. 



3 DLA-LBG CROSS-CORRELATION 

The probability of finding an object 1 in volume SVi at a 
separation r from a randoml y chosen objec t 2 can be written 
as 5P = m [1 Ci2(t-)1 5Vi (|Peebleslll98"oh . The joint prob- 
ability of finding an object 1 in the volume 1 (5Vi) and an 
object 2 in the volume 2 {5V2) at a separation r is defined 
as 5P = 711712 [1 -f (,i2{r)]5V\5V2, where 711 and 712 are the 
mean number densities of the two population. For the cross- 
correlation, we replace object 1 and 2 with DLA and LBG, 
then the joint probability between DLA and LBG is 



5P = 7?,DLA7lLBG [1 + CdLA-LBg('") (»")] 5Vdla5VlBG, 



(1) 



where tidla and 71lbg are the mean number densities of 
DLAs and LBGs, and ^dla-lbg('") is the cross-correlation 
function (CCF). 

T o estimate £dla-lbg(''), we use the method of 
iLandv fc Szalavl (|l993l ). whose variance is effectively Pois- 
son: 

Cdla-lbg(7') ~ 
-Ddla^lbg — Cdla^^lbg — -Rdla^^lbg + -Rdla-Rlbg 



-RDLAi?LBG 



(2) 



where -Ddla^lbg is the number of pairs between the two 
data samples of DLAs and LBGs separated by a distance 
r ± 5r, and likewise for other terms. The notation "-Rdla", 
for example, represents the DLA sample that has random 
coordinate positions but with equivalent number density as 
the original data sample "Ddla" . 

The method of id entifying DLAs i n our s imulations is 
described in detail in INagamine et al] (|2004al ). Briefly, we 
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Figure 1. DLA-LBG CCFs at 2 = 3 calculated with the regular unweighted method (Equation [2J . The variance of CCFs computed with 
10 different random seeds is shown with the cyan shade, and the open square symbols are the mean of 10 trials. The blue dashed line is 
the least-squa re fit to the open square points. The red solid line and the short and long dashed lines are the angular and 3-D best-fitting 
power-laws of lCooke et al.l l(2006a .b.). respectively, and the yellow shade is their 1-ct error range for the angular CCF. 



set up a cubic grid that completely covers each dark mat- 
ter halo, with the grid-cell size equivalent to the gravita- 
tional softening length 'e' of each run. We then calculate 
the Hi column density A^hi of each pixel (i.e., a grid-cell 
on one of the planes) by projecting the Hi mass distribu- 
tion, and identify those that exceed the DLA threshold of 
Nni > 2 X 10^° atoms cm ^ as 'DLA-pixels'. This method 
allows us to quantify the DLA cross-section 'ctdla' of each 
halo, and the number of DLA-pixel is Nf'^^ — ctdla/e'^- 
Here we focus on the correlation sig nal at r > 0. 4 hr^_yipc, 
because this is the scale probed bv lCooke et al.l (|2006al lbh. 
Therefore in this paper we are only concerned about the 
overall halo positions and not the exact locations of individ- 
ual DLA-pixel within each halo. The (ToLA-weighted CCF 
will be discussed in Section [3] 

First, we select the LBGs that are brighter than 
i?AB=25.5 magnitude in the D5 a nd G5 runs. There are 
30 (4030) LBGs in the D5 (G5) run. iNagamine et al.1 (|2004| ) 
have shown that the brightest galaxies with 7?ab < 25.5 in 
our simulati ons satisfy the U„G R color selection criteria for 
LBGs (e.g.. ISteidel et al.lll999l ). There are 22616 (25683) 
DLA halos with (Jdla > in the simulated volumes of 
the D5 (G5) run. The 'random' catalogues of LBGs and 
DLA halos with random positions were created with a ran- 
dom number generator from Numerical Recipes (jPress et al.l 
1 19921 ). The selected LBGs were paired with DLA halos, 
and the number of pairs that reside in spherical shells 
of [log r, log r + A log r] were counted. The maximum pair 
separations probed for the D5 and G5 runs are 10 and 
35 Mpc, respectively, with 20 bins in a logarithmic scale 
of distance r. The periodic boundary condition was taken 
into account appropriately, and the pair-search was extended 
to the next adjacent box where needed. 





Unwoif 


jhted 


0"DLA 


weighted 


Run 


To 


7 


To 


7 


D5 


2.68 


1.46 


3.35 


1.70 


G5 


2.92 


1.68 


3.30 


1.69 



Table 2. Best-fitting power-law parameters of unweighted and 
(TQLA.-weighted DLA-LBG CCFs at 2 = 3. The correlation length 
ro is in units of Mpc. For comparison, Cooke (private com- 
munication) obtained ro = 2.91 ± 1.0 Mpc and 7 = 1.21^^ ^ 
for the 3-D CCF calculated with spherical shells, and lCooke et all 
ll2006bl) reported ro = 3.32 ± 1.25 h"^ Mpc and 7 = 1.74 ± 0.36 
for the angular CCF. 



Figure [T] shows the DLA-LBG CCF computed with 
Eq. ((2]). We perform a least-square fit of the measured val- 
ues with a power-law ^(r) — {ro/r)~' , and find best-fitting 
parameters equal to (ro Mpc], 7) — (2.68,1.46) and 
(2.92, 1.68) for the D5 and G5 runs, respectively. The fits are 
shown by the blue long-dashed lines (see also Table [2]). We 
repeat the calculation of the CCF ten times using different 
seeds for generating the random positions for the 'random' 
sample to examine the statistical variance of the measured 
CCF. The variance of 10 trials is shown as the cyan shade, 
and the average of ten trials is shown with the open square 
data points. The red solid line and the yellow shade represent 
the best-fitting result (ro = 3.32± 1 25 an d 7 = 1.74± 0.36) 
and the 1-a errors of ICooke et all (|2006al lbl) from their an- 
gular CCF result. The result of the G5 run agrees well with 
that of Cooke et al.'s, and its variance is small owing to a 
larger sample than in the D5 run. The result of the D5 run 
is somewhat shallower than that of the G5 run, which could 
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simply owe to relatively small sample of LBGs in D5 and its 

sma ll box-size. 

ICooke et al.l (|2006al lbh published only the angular 
CCFs. However, they can also estimate the 3-D radial CCF 
using redshift information. The best-fitting parameters to 
the radial CCF by Cooke (private communication) using 
spherical shells is To = 3.39±1.2 Mpc and 7 = 1.61±0.3, 
which is shallower than the angular CCF r esults. As we will 
furthe r discuss in Section [S] the method of lAdelberger et al.l 
adopts cylindrical shells at small distances, which 
have larger volumes than spherical shells. This allows the 
cylindrical shell method to contain a larger number of LBGs 
an d DLAs at small dis tances, yielding a slightly steeper 7 
m ICooke et all l|2006bl ) compared to the above spherical 
shell case. We regard the comparison to the angular CCF of 
ICooke et al] (l2006bh as the primary one, because Cooke et 
al. a rgue that the angular CCF calculated by the method of 
lAdel bcracr ot al. (2003) is more robust than the 3-D radial 
calculation with spherical shells, and the values of (to, 7) 
derived from both CCFs should be equivalent theoretically 
(see Section [S}. 



4 (TDLA-WEIGHTED CCF 

In Section [31 we calculated the CCF assuming that there is 
one DLA per halo. This assumption is valid as long as we 
are concerned with the CCF at scales of r > 300/i~^kpc, 
because in the real observations the same dark matter 
halo is rarely probed b y a close pair of quasars. However, 
iNagamine et al.l (|2004bl . Fig. 1) showed that the DLA clouds 
have extended distributions in massive dark matter halos. 
Therefore, it may be more desirable to take the DLA cross- 
section of each halo into account when calculating the CCF. 
A simple way to achieve this is to weight the number of DLA- 
LBG pairs by the number of DLA-pixels of each halo. Since 
the displacement between DLA-pixels within a single halo is 
typically much smaller than the distance between LBG-DLA 
pairs, we do not count the individual pairs between LBG 
and DLA-pixels. Instead, we treat it as if all DLA-pixels 
are located at the halo center, and weight each DLA-LBG 
pair-count by the number of DLA-pixels Ni (hereafter we 
drop the superscript 'DLA' for simplicity) and compute the 
o"DLA-weighted CCF as 

Cdla-lbg('') = 

iVi-DDLADLBG-iVjDDLA-RLBG-^'i-RDLADLBG+^fi-RDLA-RLBG 



JVi-RDLA«LBG 



(3) 



For the 'random' DLA dataset, we shuffle the original Ni 
list randomly and make new pairs with different DLA halos. 
Again, 10 realisations of the random dataset have been used 
to examine the statistical variance of the estimated CCF. 

The results for the (jDLA-weighted CCF is shown 
in Figure (2] We find best-fitting parameters equal to 
(ro [/i"^Mpc],7) = (3.35,1.70) and (3.30,1.69) for the D5 
and G5 runs, respectively, as shown by the blue long-dashed 
line (see also Table [2]). Both results show remarkable agree- 
ment with the best-fitting values of iCooke et al.l (|2006bl . 
ro = 3.32 ± 1.25 and 7 = 1.74 ± 0.36). The result of D5 
is somewhat noisy at r < 1 Mpc, which originates from 
the noisy pair-count of A^j-DdlaOlbg- 



The parameter values given in Table|2]clearly show that 
the (TDLA-weighted method gives larger values of ro and a 
slightly steeper power-law slope. In a CDM universe, the 
number of low-mass halos is far greater than that of mas- 
sive halos. Therefore, even a small weighting by Ni boosts 
up the overall pair-count, yielding a stronger correlation sig- 
nal compared to the unweighted case. The larger LBG sam- 
ple in the G5 run makes its result more robust against the 
weighting procedure than that of the D5 run, and the slope 
7 in the G5 run does not change very much between the two 
calculation methods. 



5 ANGULAR CROSS CORRELATION 
FUNCTION 

In observational studies, a different method is usually used 
to obtain the values of (ro, 7) compared with what we de- 
scribed in Sections [3] and |4j because the precise estimation of 
any LBG position along the line of sight is difficult to achieve 
owing to redshift uncertainties caused by peculiar velocities 
and galactic winds. With such imprecision, it is not possi- 
ble to measure the CCF at scales r < lh~^ Mpc reliably. 
Therefore, rather than attempting to estimate the 3-D dis- 
tance between DLAs and LBGs, observers usually employ 
the angu lar CCF us ing the projected data on the sky. For 
example. ICooke et a l. ( 2006aJ3) computed the an g ular C CF 
using the method proposed by lAdelberger et all l|2003h . In 
order to compare our results with those by Cooke et al's, we 
briefly describe the calculation method of lAdelberger et al.l 
I2OO3) , and then describe how we perform our measurement 
of the angular CCF. 

With a power-law assumption, the expected number of 
pairs for the projected angular CCF is 



i^pije < r^) = 



2r;, 



(4) 



where B a nd I:c are the beta and incomplete beta functions 
with (;e.g.. |Press et al.lll992l ) 

a; = r^ (r^ -l-rfl) \ (5) 

lAdelberger et al.l (|2003l ) proposed to count the number of 
pairs in cylindrical shells of angular separation re ± Sre 
and redshift separation r^ ± Sr^, rather than using spher- 
ical shells. By setting to 



r, = max | 1000 kms"^ ^i^r^, 7re 



H{2 



(6) 



the lower limit ensures that the redshift errors do not lead 
to the underestimate of the number of pairs, and the upper 
limit allows sufficient distanc es to include enough correlated 
pairs l|Adelberger et al.ll2003l ). 

For our calculations, we focus a,t z — 3 and thus r^ — 
max(12.8 /i^^ Mpc, 7re). With simple algebra. Equation Q 
can be converted to a more familiar power-law form: 



are) 



2r„ 



re 



B 



i 1. 

2' 



i 1. 

2' 



(7) 
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Figure 2. DL A-LBG CCF at z = 3 calculated by the troLA-weighted method (Equation [HJl. The yellow shade shows the upper and 
lower limits of lCooke et al. I ll2006allbh 's best-fitting power-laws. The variance of the CCF using 10 random seeds is shown with the cyan 
shade. The blue dashed lines are the best-fits for this work, and th e red solid line and the short and long dashed lines are, respectively, 
the angular and 3-D best-fitting power- laws of lCooke etHI (l2006al lbl) . 



where is set to r-max- We change from spherical coor- 
dinates to cylindrical coordinates, and set the number of 
cylindrical bins to 20 in a logarithmic scale as before. All 
pair searches are extended to the adjacent box using peri- 
odic boundary conditions, if appropriate. 

A few assumptions must be made while we deal with 
the beta and incomplete beta functions. There are two pa- 
rameters (7 and x) that must be given to calculate the val- 
ues of B and I^^. To calculate 7, we first plot Equation (O 
without B and Ix (i.e., 2rraaxi^p{rg)/rg) and find the best- 
fitting value of 7. The value of x is determined by Vz and 
re as shown in Equation ((5|. By setting ~ ^max, the an- 
gular separation will be divided into two different regimes. 
Within the smaller angular separation range (100 kpc < 
rg < 1.83 Mpc), the correlated pairs are counted up to 
the maximum radial distance of rmax ~ ±12.8 Mpc for 
a cylinder centered on an LBG or DLA, while in the larger 
separation range [rg > 1.83 Mpc) all the correlated pairs 
within Tniax = ±7re are counted. With the fixed value of 7 
obtained above and 20 different values of x for each bin, B 
and Ix can be calculated for each bin. 



The angular CCF results of our calculations are shown 
in Figures |3] and |4] for both the unweighted and the ctdla- 
weighted method. The best-fitting power-law parameters are 
given in Tab l e El A gain, the agreement with the results of 
ICooke et al.l (|2006al lbh is quite good. Similarly to the 3-D 
CCF case, the crDLA-weighted case gives a slightly larger ro 
and steeper 7 than the unweighted case. The unweighted 
case of D5 is shallow with 7 — 1.42, but in the ctdla- 
weighted case, 7 ~ 1.7 is recovered. 



Unweighted iTDLA-weighted 



Run 


To 


7 


ro 


7 


D5 


3.44 


1.42 


3.94 


1.71 


G5 


3.33 


1.69 


3.75 


1.70 



Table 3. Best-fitting power-law parameters for the angular CCF 
at 2 = 3. The units of the paramet e rs are the same as in 
Table [2] For comparison, ICooke et all l l2006bl) reported To = 
3.32 ± 1.25 h-^ Mpc and 7 = 1.74 ± 0.36 for their angular CCF. 

6 AUTO-CORRELATION FUNCTIONS 
6.1 LBG auto-correlation 

The auto-correlation function (ACF) also gives important 
constraints on the distribution of the population under 
study. In this section, we calculate the 3-D LBG ACF by 
changing all subscripts in Equation ((2)1 to 'LBG': 

^lbg-lbg('") = 

I^LBG^^LBG — 2DlBG-RlBG + i^LBG^^LBG 
RLBcRhBG 

(8) 

Our result for the ACF is shown in Figure (5] and the 
best-fitting power-law parameters are ro = 3.98 Mpc 
and 7 — 1.55. The last two data points were not included 
for the power-law fit because they are likely underestimated 
owing to the limited box-size. Our values of rp and 7 agree 
well with the obse rvational e stima tes of lAdelberg er et all 
(2003) and Ad elberger et al.l (120051 ). who measured values 
of To = 4.0 ± 0.6 h-^ Mpc and 7 = 1.57 ± 0.14 for the LBG 
ACF at z ~ 3, with a correction for the integral constraint. 
This correction owes to the finite size of the observed field- 
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Figure 3. DLA-LBG angular CCFs at 2 = 3 computed by the unweighted method for the D5 and G5 runs. Other features are the same 
as described in Figure [T] 
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Figure 4. DLA-LBG angular CCFs at 2; = 3 computed by the CDLA-weighted method for the D5 and G5 runs. Other features are the 
same as described in Figure [Tl 



of-view, and it must be added to the computed correlation 
function as follows: 



ACF value. Therefore, we neglect the correction by IC in 
this paper. 



C'(r)=CM+IC, 



(9) 



where ^'(r) and ^(r) are the corrected and computed ACF, 
respectively, and IC is the integral con s traint . Foll owing the 
metho d described in lAdelberger et al] l|2005h and lLee et al.l 
(|2006l ) ■ we calculate the value of IC and find that it is not 
significant in our simulations compared to our computed 



The dark matter ACF (the red filled triangles in Fig- 
ure [3 was also computed as described in iNagamine et alj 
(120081 ) in order to calculate the bias of LBGs against the 
dark matter distribution (see Section [7]). 
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Figure 5. LBG auto-correlation function at z = 3 for the G5 run. 
Tlie y ellow shade shows tiie l -cr range of the best-fitting power- 
law of lAdelbereer et al.l |200i). The variance of the ACF using 10 
random seeds is shown with the cyan shade. The red solid and blue 
dashed lines are the best-fitting power-laws of lAdelberger et al.l 
1I2OO5I ) and this work, respectively. The last two data points were 
not included for the power-law fit. The red filled triangles show 
the dark matter ACF at the same redshift. 



6.2 DLA auto-correlation 

Similarly to the LBG ACF, it would be useful to compute 
the DLA ACF in order to estimate the DLA host halo mass. 
In this section, we calculate the DLA ACF with both the 
unweighted and the (JDLA-weighted methods. By replacing 
all subscripts to 'DLA' in Equations ((S)) and © , we obtain 



CDLA-DLA(r') = 

OdlaI^dla 



2DdLA^JdLA + ^?DLAi?DL/ 



(10) 



and 



AWoightcd f \ 
SDLA-DLAV ) 



^DLA^DLA" 



^DLA-^DLA 



+ NiNjIi 



DLA^DLA 



(11) 

where NiNjD}yi^j^D^j^j^ and NiNj D^^j^p^R^j^j^ are the num- 
bers of data-data pairs and data-random pairs, weighted by 
the number of DLA pixels Ni and Nj . As before, ten different 
realisations of random dataset have been used to examine 
the statistical variance. 

Our DLA ACF result is shown in Figure [6l and 
we find the best-fitting power-law parameters of To = 
2.43/i"^Mpc and 7 = 1.60 for the unweighted ACF, and 
To = 2.99/i~^Mpc and 7 — 1.55 for the (TDLA-weighted 
ACF, as summarized in Table |4l The values of 7 are similar 
to those for the LBG ACF with 7 ~ 1.6, but To is much 
smaller. This is owing to the lower average DLA halo mass 





To 


7 


LBG-auto 


3.98 


1.55 


DLA-auto (unweighted) 


2.43 


1.60 


DLA-auto ((Tj3LA-wcighted) 


2.99 


1.55 



Table 4. ACFs of LBGs and DLAs for the G5 run. The results of 
unweighted and crnLA-wsighted methods are g iven for the DLA 
ACF. To is in units of Mpc. For comparison. lAdelbereer et al.l 
l|2005l ) reported ro = 4.0 ± 0.6/1"^ Mpc and 7 = 1.57 ± 0.14 for 
the LBGs at z ~ 3. 



compared to the LBG host halos, as we will discuss further 
in Section [T] 



7 BIAS AND HALO MASSES 

Comparing the correlation functions of DLAs and LBGs 
with that of dark matter gives the measure of 'bias' for the 
spatial distribution of these populations against that of dark 
matter. Figure [7] shows the bias, defined as 6 = \/^77Cdm, 
as a function of distance r, where i = LBG or DLA. This 
definition is based on the linear bias model. 



(12) 



The corresponding expression for the cross- correlation is 
(|Gawiser et al.ll200ll ) 

^DLA-LBg('") = feDLA6LBG?DM('")- (13) 

Therefore, the two lines for the CCF in Figure [7] are in fact 
showing \/6dla6lbg, as indicated on the axis on the right- 
hand- side. Taking the ra tio of the above two expressions 
gives (|Cooke et al.ll2006bh 



CDLA-LBG(r) _ feDLA 



(14) 



ClBG &LBG 

In all cases shown in Figure [T] the bias slowly decreases 
with increasing distance. The upturn at r = 20/i~^Mpc 
for the LBG ACF is probably just noise. We take a sim- 
ple average of bias values across the logarithmic bins at 
r = 1.40 - 14.5 hr^ Mpc, and find 6 = 2.60, 2.47, 2.30, 2.24 
and 1.93 for LBG ACF, DLA-LBG CCF (ctdla- weighted), 
DLA ACF (cTDLA-weighted), DLA-LBG CCF (unweighted), 
and DLA ACF (unweighted), respectively. The values of 
To also reflect the sizes of average bias values. We took 
the above range of scales for taking the average because 
most of the recent observations are probing the scale of 
r ~ 1 - 10/i~^Mp c. 

Gawiser et al.l (|2007f) used the results of 

lAdelberger et al.l (|2005h to obtain an average bias of 
&LBG = 2.5 ± 0.4 for LBGs at 2; ~ 3. Our average bias 
va lue of 2.60 for the LB G ACF is very close to that 
of lAdelberger et all (|2005h . and at the lower end of the 



estimate of 6lbg = 3.0 ± 0.5 bvlLee et all (l2006h 



The model of ISheth fc TormenI (|l999l ') shows that an 



understanding of the unconditional mass function can pro- 
vide an accurate estimation of the large-scale bias fac- 
tor. From our average bias, we calculate the mean halo 
mass for LBGs and DLAs (using the unweighted and the 
cpLA-weighte d results) based on the method described in 
IMo fc Whit3 12002,), as shown in Tabled Our calculation 
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Figure 6. DLA auto-corrclation function calculated with unweighted and aj^i^^-weighted metliod for tlie G5 run. The variance of ACFs 
using 10 random seeds is shown with cyan shade. The blue dashed lines are best-fits for this work. 



of LB G halo mass is very close to that by I Adelberger et al.l 
HqoI), MiJ;^^ = 10"-^ - 1O^^-*M0 (yellow shade inFig. 0, 
which is very encouraging. Finally, iBouche et all l|2005t l es- 
timated (logMnLA) = 11.13 ± 0.13 from observations and 
(logMoLA) = 11.16 from simulations. These values are 
somewhat higher than the upper limit of our unweighted 
DLA halo mass an d very close to our croLA-weighted one. 
ICooke et al.l ()2006ah also obtained a similar value of Mhaio — 
1O"-^M0. 

Alternatively, we can directly calculate the mean 
DLA halo mass using the simulation result without going 
through the bias argument. For the G5 run, the mean is 
^og{M^tio) = 11-5 and (logA^^^ ^ ^hese val- 

ues are somewhat higher than the mean halo mass re- 
ported in Table [5] However, the values of (Mhaio) in Ta- 
ble [5] are computed from the average bias within the range 
of r = 1.40— 14.5 Mpc, and they could become higher if 
we included the bins at smaller scales. Since observers probe 
mostly r ~ 1 — 10 Mpc, the values reported in Table [5] 
are more appropriate for co mpari son with observations. 

iBouche &: Lowenthall ()2004l ) defined the parame- 
ter a as the ratio of correlation functions: a = 
6ccf(Mdla)/&acf(A^lbg)- If the value of a is larger (or 
smaller) than unity, then the mean halo mass of DLAs is 
more (or less) massive than that of the LBGs. The ratio 
of the average bias of LBG ACF and DLA-LBG CCF is 
Q = 0.742 for our results. This value is in good agree- 
ment with the observationa l estimates of a — 1.62 ± 1.32 
jBouche fc Lowenthajp003). a = 73 ± 0.08 (|Bouche et all 
l2005l ). and a ^ 0.771 (|Mo fc Whitd[200^ ). 





bias 


log(Afhalo> 


LBG-auto 
DLA-auto (unweighted) 
DLA-auto (tTDLA-WEighted) 


Z.DU_Q 22 

^■''•-'-0.13 
o qri-t-0.23 


ii.48«:^2^ 

iU. (U_Q j^g 

11 16+°-23 

il.iD_Q j^y 



8 DISCUSSION AND CONCLUSIONS 

Our study represents a first attempt to calculate the DLA- 
LBG cross-correlation function at z; = 3 using cosmologi- 



Table 5. Average biases and halo masses of LBGs and DLAs 
for the G5 run. The plus and minus values next to the average 
bias show the upper and lower limits at 1.40 < r < 14.5 Mpc. 
Mean halo masses a re computed from the second column following 
IMo fc Whit3 1 I2OO2I) and given in units of Mq. 



cal SPH simulations. We calculated the DLA-LBG CCFs in 
several different approaches: 3-D, angular, unweighted, and 
o"DLA-weighted. We also computed the auto-CF of LBGs and 
DLAs, and the bias against dark matter. In c omparison to 
observa tional data bv I Adelberger et al.l (|2005t ): [Cooke et al.l 
(|2006al lbl). we find good agreement between our simulations 
and observational measurements. Our results suggest that 
the spatial distribution of DLAs and LBGs are strongly cor- 
related. 

Let us summarize some of the main conclusions of this 
work. In the first part of this paper, our results on the 3-D 
CCF calculated with spherical shells (Table[21) are to be com- 
pared with the 3-D spherical shell result by Cooke (private 
communication), ro — 3.39±1.2 Mpc and 7 — 1.61±0.3. 
Our results are consistent with Cooke's within the error. The 
shallow slope of Cooke's above estimate probably owes to the 
limited sample size in the spherical shell at small distances, 
as we discussed in Sections [3] and [S] 

In the second part, we have replaced the spheri- 
cal shell method with the projected approach used in 
lAdelberger et al.l (120031 ) and lCooke et all l|2006br ). and cal- 
culate the best-fitting values given in Table [S] Encourag- 
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Figure 7. The biases of all correlation functions at z = 3 that 
we computed in this paper for the G5 run. The tick marks on 
the left-hand-side s how the host halo m asses calculated with the 
method described in IMo fc Whitj 1I2OO2I). The yellow s hade shows 
the upper and lower limits bv" ^^e^erger et'aJ ] l|2005l 'l. 



ingly, our results are within the u pper and low er limits of 
the observational measurement bv lCooke et al. (2006a b). 

Finally, we also analyzed the auto-correlation functions 
of LBGs and DLAs at « = 3 (Table |4| found in our sim- 
ulations. Our results for t he best-fittinR pa r amete rs of the 
LBG ACF agree well with I Adelberger et all lj2005l 'l. We cal- 
culated the integral constraint and confirmed that it was not 
significant. Our result show that LBGs are more strongly 
correlated than DLAs, and have higher mean halo mass. 

Figure [8] summarizes the best-fitting power-law param- 
eters for all the correlation functions that we obtained in 
the earlier sections. In most cases, the slope 7 falls into the 
range « 1.5 — 1.7 and the variation is not very large. The 
correlation length To shows a much larger variation from 
2.5 ft"'^ Mpc to 4/i~^Mpc, depending on the sample and 
calculation method. This trend is similar to that seen by 
ICooke et al.l l|2006bl . Fig. 8). In general, the (TDLA-weighted 
method gives a larger ro than the unweighted method, and 
the D5 run tends to underestimate the value of 7 owing to 
the small LBG sample size. 

Finally, the LBG bias, derived from the LBG ACF in 
Section [3 has led to the upper and lower limits of the 
LBG dark matter halo mass of log{Afhaio) = 11.48lg'22 (s^e 
Table [5]) . This result is consistent with observational esti- 
mates of the LBG ha l o mass of Mhaio ^ IO^'^Mq, (e.g., 
lAdelberger et al.l 1 19981 : ISteidel et all 119981) and within th e 
limit of Mhaio = 10"'^ -10"-^Mm (lAdelberger et al.ll2005h . 
Similarly, we derived the DLA biases, and obtained the me an 
DLA halo masses as shown in TablelSl lCooke et al.l (|2006al )'s 
measurement showed a DLA galaxy bias of ^dla ~ 2.4 and 
an average DLA halo mass of Mhaio ~ 10^^'"^ Mq. Our av- 
erage DLA bias and halo mass estimates are in good agree- 
ment with theirs. We also examined the ratio of bias values 



Figure 8. Summary of best-fitting power-law parameters for all 
correlation functions that we obtained in earlier sections. 



defined as a = 6ccf /6acf (|Bouche fc Lowenthalll20o3 ). and 
found that our value of a = 0.742 agrees well with the ob- 
servational estimates. This again shows that the mean halo 
mass of DLAs is less than that of the LBGs. 



The fact that 



(AfLifG) is greater than {M^^^t) 



a natural outcome because the LBG sample is limited 
to the bright star-forming galaxies with Rab < 25 and 
M* ~ 10"' - 10" Mq, whereas the DLA Hi gas is 
present in numerous lower mass halos below the LBG 
threshold. Given the good agreement between our re- 
sults and the observations, and considering all the ac- 
cumulated ev idence that sugg e sts a high halo m a ss for 
LBGs (e.g., lAdelberger et al'f Il998l: iBaugh et all Il998l: 



iGiavalisco et al.' ' 



1998 



Katz et al. ' '1999'r Kauffman n et al 
1999; Papovich et al 



20031) no 



I2OOII : [Shaplev et al.1 l200il : ISteidel et all 1 19981 ) . we con- 
clude that scenarios where the majority of LBGs are 
merger-induc ed starburst systems associated with low- 
mass halos jLowcnthal ot al.' 1997; 'Sawicki fc Ye3 Il998l : 
ISomerville et al. 2001; Wcatliprlcy &_ Warren, 
longer appear to be viable models for LBGs. 
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